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Abstract 

In this paper, we consider the problem of Multiuser 
Detection (MUD) in Direct Sequence/Code Division Multiple 
Access (DS/CDMA) communication systems. The optimum 
detector for MUD is the Maximum Likelihood (ML) detector, 
but its complexity is very high and involves an exhaustive 
search to reach the best fitness of the transmitted and 
received data. Thus, there has been considerable interest in 
suboptimal multiuser detectors with less complexity and 
reasonable performance. If we look to the problem from an 
optimization viewpoint, what makes our work difficult is 
high interaction and Epistasis among variables which should 
be optimized. We propose to work on the problem and try to 
reduce its complexity instead of using more complex and 
powerful optimization algorithms. To do so, we analyze 
problem carefully and propose a new method to overcome 
Epistasis. In this paper, by using sign detector and longer 
code lengths in DS/CDMA system, a simple detector will be 
proposed. This method is compared with ML and two 
model-free optimization methods: Particle Swarm 
Optimization (PSO) and Ant Colony Optimization (ACO) 
algorithms which have been used for MUD in DS/CDMA. 
Simulation results show our method dramatic efficiency. 
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Introduction 

In a DS/CDMA system, the receiver is a Matched Filter 
Bank (MFB), which comprises the conventional 
detector (sign detector). This type of receiver is unable 
to optimally recover the signal, when the channel is 
contaminated by additive white Gaussian noise 
(AWGN), and suffers from flat or frequency selective 
fading; because the DS/CDMA signal is affected by 
Multiple Access Interference (MAI) and also by the 
Near-Far Ratio (NFR) [1]. In fact, the signature signals 



of different users are not completely orthogonal to 
each other, and cross correlation among these signals 
results in multiple access interference. Therefore, the 
conventional matched filter detector, as in single user 
communication, is no longer effective and causes 
many problems [2]. In 1986, Verdu in [3] proposed the 
optimum multiuser detector (OMUD) which consists 
of a bank of matched filters followed by a maximum 
likelihood sequence estimator (MLSE). The MLSE 

detector generates a maximum likelihood sequence, ^ , 
which is associated with the transmitted sequence, as 

presented in Fig. 1 [1]. The vector ^ is estimated in 

order to maximize the sequence transmission 

r(t) r(t) 
probability given that v J is received; where v J is 

extended for all messages, considering all the 

transmitted messages with the same transmission 

probability [1]. The OMUD has a computational 

complexity which grows exponentially with the 

number of users. Thus, since CDMA systems could 

potentially have a large number of users, the OMUD is 

impractical to implement for them. Therefore, many 

researches have focused on sub-optimum detectors 

with less complexity and a performance which is 

almost as high as the OMUD. Alternative detectors for 

OMUD include the Decorrelator proposed by Verdu in 

[3] and the MMSE detector recommended by Poor and 

Verdu in [4]. These algorithms have reasonable 

computational complexity, and their performance is 

comparable to that of the optimum receiver, but they 

yield a degraded communication system in the sense 

ofBER [5]. 

According to the problem of ML detector, many 
methods for suboptimal detection have been proposed. 
Some heuristic methods have been developed, such as 
genetic algorithm (GA). The first GA-based multiuser 
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detector (GA-MUD) was proposed by Juntti et al. [6] 
which assumed a synchronous CDMA system model 
communicating over an AWGN channel. After that in 
[7, 8] the new approach of GA is proposed. Tabu 
search algorithm [9] and simulated annealing 
algorithm (SAA) [10, 11] are the other new approaches 
of Multiuser detectors. 

In fact, all the approaches for MUD are classified into 
two global approaches: the first is the proposition of a 
new suboptimal Multiuser detector with tradeoff 
between performance and complexity such as MUD 
based on Particle Swarm Optimization (PSO) and 
wavelet transform [12-17] and Ant Colony-based 
MUD in [18-21]. The second approach is the reducing 
of MAI and applying simple conventional detector. In 
other words, instead of using more complex methods 
for multiuser detection, we try to reduce Multiple 
Access Interference in DS/CDMA system and then 
apply simple detector such as conventional detector. 
For this purpose, we must introduce Epistasis. 

Epistasis is a phenomenon where the effects of one 
gene are influenced by one or several other genes. The 
word "Epistasis" was used by William Batson [22] in 
genetic science to describe a gene that affects the 
operation of other genes. 




FIG. 1 BASEBAND DS/CDMA BLOCK DIAGRAM, RECEIVER 
WITH MUD [1]. 

In optimization literature, Epistasis is used with 
similar definition: the interaction between 
optimization variables [23]. Highly Epistatic problems 
are very more difficult to be optimized. In fact in every 
nontrivial optimization problem, there should be at 
least a little degree of Epistasis; otherwise, the problem 
will be trivial [24]. Since in such cases, the N-variable 
non-Epistatic function can be divided into N separate 
functions each with one variable. In our optimization 
problem, MAI is the cause of Epistasis. So to design a 



suitable detector for DS/CDMA, we try to reduce 
Epistasis (here, MAI) and use conventional detector 
rather than using more complex Multiuser detectors. 

MAI is the result of correlation between signature 
signals of users. Correlation could be shown in square 
matrix form in which elements on main diagonal are 
autocorrelations and elements outside main diagonal 
are cross correlations between the signature signals of 
users. For M-sequence codes if code with higher 
length is used, it will be claimed that autocorrelation 
elements will remain constant but cross correlation 
elements decrease proportional to inverse of code 
length [25], i.e. longer M-sequence codes, less MAI. It 
is important to note that long codes have some 
advantages such as interference randomization and 
improving the system capacity. Based on the above, 
our analysis focuses on DS/CDMA systems with long 
spreading codes. Our proposed method is compared 
to the ML detector, Matched Filter detector, Ant 
Colony and PSO-based detectors, and we show that 
the performance of the proposed method is better with 
less complexity. 

The remainder of this paper is organized as follows: 
Section 2 describes our asynchronous CDMA system 
and section 3 highlights the simple method used to 
implement our proposed detector. The simulation 
results are presented in Section 4, while Section 5 
provides comparison of the complexities associated 
with our algorithm and several effective methods. 
Finally, some conclusions are dawn in Section 6. 

System Model 

In a DS/CDMA system with binary phase-shift keying 
modulation (BPSK) shared by k asynchronous users, 
as illustrated in Fig. 1, the k-th user transmitted signal 
is given by [1]: 

x k (0 = V^Z^\ {t-iT b )cos(a) c t). 

(1) 

Where Pk represents the k-th user transmitted power; 

^(O T co 

k is the i-th BPSK symbol with period b ; c is the 

s (t) 

carrier frequency and k w corresponds to the 

[or ) 

spreading sequence defined in the interval L ' b ' : 

**(0=ZW- nr <:)W 0<t<T„ 

n=0 ( 2 ) 
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Where 



c*,„e{-Ul} 



is the n-th chip of the sequence 
T 

with length N used by the k-th user; c is the chip 

T 

period and the spread spectrum processing gain, c is 

equal to N ; the pulse shaping is assumed 

rectangular with unitary amplitude in the interval 

[or ) 

L ' c ' and zero outside. 

Assuming a frame with / bits for each user 
propagating over L independent slow Rayleigh 
fading paths, the baseband received signal in the base 
station is [1]: 



I -l K L 



i =0 k =1 / =1 

(3) 

Where K is the number of active users, t G \P^ b ) . 
The amplitude k is assumed constant for all / 

is the transmitted 



transmitted bits, ^ E t ^' + ^ 



information bit, k denotes a copy of the signature 

sequence assigned to the k-th user and Tk represents 
the random delay associated to the k-th user; the 
complex low-pass impulse response of the channel for 
the k-th user over the i-th bit interval can be written as 
[1]: 

(4 ) 

X a {i) 
Where kJ is the propagation delay and kJ is the 

complex channel coefficient whose amplitude has 

Rayleigh distribution, and its phase is uniformly 

distributed over ^' . Finally, w(t) represents the 

AWGN with bilateral power density equal to 2 . 

Using vectorial notation, Equation. (3) can be stated as 
[1]: 



r{t) = ^ T {t-iT b )Aa {i) b {i) +w (t). 



(5) 



A =diag\AJ,AJ 9 ...,A t l] . , 
Where 6 L 1 2 * J is the diagonal 

matrix for the users' amplitude including the path 

losses and shadowing effects, and ^ LxL is the identity 
matrix with a dimension equal to L ; s is the vector 
of users signature sequence, and a is the diagonal 
channel gain matrix as: 

«<'>=^[ fl a>,..,<\aa,..,a«>,..,fl«i] 



And the data vector is given by: 



,(0 - 



b (i) b (i) b (i) 
u \ -> u 2 ->""> U K 



(6) 



(7) 



representing the 1><L k-th user bit vector. 

If we use the conventional Rake receiver which 
consists of a bank of KL filters matched to the users' 
signature sequence [2], then the output for the i-th bit 
interval can be expressed as [1]: 



% =^r (f> t (t -iT b -r k) )dt =A k T b p$bP +/£> +4i . 

(8) 



Where the first term corresponds to the desired signal, 
the second refers to the MAI over the Lr -th multipath 
component of the k-th user, and the last represents the 
filtered AWGN. It should be noted that the auto- 
interference term is neglected for simplicity. MAI is 
the result of cross correlation between signature 
signals of users as expressed [1]: 



r j* ( Tj )=fo s j ('K (' +iT b +T ) dt 



(9) 



If we use the Rake receiver, we need to estimate some 

(t) 

parameters such as channel coefficients and delay v ' . 
When the number of users increases, the significance 
of interference rises and owing to this fact, the 
conventional detector performance in MUD is 
degraded. As mentioned in the previous section, we 
use the optimum detector to solve this problem. The 
optimum detector is the maximum likelihood 
sequence detector that selects the most likely sequence 
of transmitted bits given by the observations at the 
receiver. In this context, the K-user, L-paths, I-frame 
and asynchronous channel scheme can be viewed as a 
KLI-user synchronous channel scheme, and then the 
KLI-user vector B can be written as [1]: 
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B = 



_ b 



(I A)' 



(10) 



Based on [27], Verdu proved that in order to select the 
maximum likelihood sequence B, we must maximize 
the log likelihood function (LLF): 



/ (b) = 2Rq{b t a H Ay}-B J aARAa H B . 



(11) 



Where R is the cross correlation matrix and y is the 
output vector as [1]: 



y 



y (0) ,y m 



(12) 



Neglecting the channel effect, we can state (11) in a 
simple form as [1]: 



/ (B) = B T r-B T RB . 



(13) 



Where r is the received signal and B is the 
transmitted sequence to be guessed and we neglect 
other parameters in (11) due to the channel effect. The 
complete frame with the estimated transmitted bits for 
all K users can be obtained through optimization of 
(13), resulting [1]: 



B = ars 



max 



[/<«)]] 



Be{-IM}" 



(14) 



The OMUD attempts to find the best vector of data 
bits but because of high complexity and unfeasible 
implementation, it is an inefficient method for 
multiuser detection. Because of the high 
dimensionality of OMUD, all suboptimal algorithms 
try to find a solution following an objective function 
which is able to improve the performance of multiuser 
detection. These attempts try to reduce the complexity 
of OMUD and maximize the DS/CDMA mean 
performance. Most efforts concentrate on approaching 
the performance of ML algorithm with less complexity 
and reliable applicability along with the least possible 
error. In the next section, we propose our algorithm to 
achieve this goal and compare our algorithm to other 
efficient algorithms available in the literature. 

Proposed Algorithm 

Sign Detector 

In conventional single user detection, the receiver for 
each user consists of a demodulator that correlates the 



received signal with the signature sequence of the user 
and passes the correlator output to the sign detector, 
which makes a decision based on the single correlator 
output. Thus the sign detector neglects the presence of 
the other users of the channel or equivalently assumes 
that the aggregate noise plus interference is white or 
Gaussian. Clearly if the signature sequences are 
orthogonal, the interference from other users vanishes 
and the conventional signal detector is optimum. 
Because it is not possible to design signature 
sequences for any pair of users that are orthogonal for 
all time offsets, thus the conventional detector is 
vulnerable to interference from other users [2]. 

Now we should choose between two options. Using 
more complex and powerful detectors to overcome 
MAI and give a suboptimum detector near ML 
detector, or try to eliminate Barriers mentioned above 
in sign detector method. 

Epistasis and MAI 

In genetic science, Epistasis is a phenomenon where 
the effects of one gene are influenced by one or several 
other genes. This concept exists in optimization 
problems too where the values of one optimization 
variables cannot be labeled (as good, bad and so on) 
without considering the status of other variables. For 
example consider the optimization test functions in 
tables 1 and 2. These functions are some of the most 
famous test functions which are commonly used in the 
literature [27-31]. In table 1 you can see some non- 
epistatic problems. In these problems you can 
optimize the objective function for each variable 
separately, while in the problems such as the ones in 
table 2, the important matter is the arrangement of the 
variables and their values beside each other. In such 
problems there is no superposition and the whole is 
something more than sum of the individuals. That 
extra comes from the information existing in the 
structure. 

TABLE 1 FUNCTIONS WITH ZERO DEGREE OF EPISTASIS 



Function 


Formula 




N 

n=l 


F 2 


x sin(4x ) + l.lv sm(2y ) 


F 3 


10Ar+X[*„ 2 -10cos(2,rx„)] 

n=\ 
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TABLE 2 EPISTATIC FUNCTIONS 



Function 


Formula 


F 4 


z{ioo[x n+1 -/] 2 + [i-^] 2 J 


F 5 


sin 2 ( \jx 2 + y 2 J - 0.5 
1 + O.ll x +y 1 



(15) proportional to 1/N, [25] i.e. longer M-sequence 
codes, less Epistasis. 



Next, let us take a closer look at the LLF. Considering 
Equation (13), the objective function is comprised of 

two terms: B T r and B T RB . p Q r maximizing the B T r 
part, considering the fact that each element of B can 
only take the values of either 1 or -1, it is enough for 
each element of B to be either 1 or -1 according to the 

sign of corresponding element in r . Since the B T RB 
term is a square matrix, if the non-diagonal elements 

were small, it would be practically equal to B T B m \ n 
this case, it is insignificant if the values of the elements 
of B are either 1 or -1. The main problem is that some 
non-diagonal matrix elements are not considerably 
small and cannot be assumed as negligible almost 
always. Because the matrix R is symmetric, these 

2R B T B 

elements add terms of the type ij 1 7 to the 
objective function. In this case if ij is positive, the 



B 



B 



fact that * and j bearing different sign helps 
increasing the value of the objective function, and if 

ij is negative, their similar signs help this fact. In 
other words, the problem variables are no longer 
independent and their interaction is influential. 

In the first glance it seems that we should try to find 
the optimization methods which have potential ability 
to overcome this problem to some extent as in [13, 15]. 
But what if we could decrease the Epistasis and make 
the problem trivial. If we do this we can solve the 
problem with sign detector easily. 

Using Longer M-Sequence Codes 

If there is not much limitation on band width in 
communication system, cross correlation between 
codes can be decreased using longer M-sequence 
codes. In other words the probability of user 
interference will be reduced and this means that 
elements that are not on the main diagonal of matrix R, 
using this work, are decreased noticeably but the main 
diagonal elements are fixed. 

For M-sequence codes if codes with higher length are 
used, cross correlation will be reduced according to 



1 if i = j 

*H 1 * • • 

In 



(15) 



So by using longer codes, Epistasis would be 
decreased and sign detector could be applied for 
detecting the transmitted signal with minimum error. 
This property is also verified for Gold codes or other 
codes that used in some communication systems for 
spreading signals. Therefore by using longer codes, 
the capacity of system is increased and the interference 
among users is decreased and consequently the simple 
and low complexity detectors can be applied for MUD. 

For better presentation of MAI decreasing, one 
simulation is done and the average of cross correlation 
of users in different code lengths is evaluated and 
depicted in Fig. 2. In this figure, it is showed that 
increasing the code lengths reduced the cross 
correlation mean value exponentially rather than 

linearly. In fact the relation is followed by ^ 
that y and x are presentations of cross correlation 
mean value and code lengths respectively. From 
simulation and using interpolation, the values of a and 
b are approximated as, b = 10.94 and a = 4. This fact 
verifies our method effectivity. Since by small increase 
in code-length a considerable decrease in MAI is 
obtained. 




FIG. 2 CROSS CORRELATION MEAN VALUE VERSUS 
DIFFERENT CODE LENGTHS, THE EXPONENTIAL RELATION. 

Numerical Results 

In this section, the performance of the algorithm 
described in Section 3 is compared to Matched filter 
receiver (MF), Maximum Likelihood method (ML) and 
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two optimization-based algorithms, PSO and ACO 
considering the BER as the main figure of merit. The 
convergence of each algorithm versus optimization 
parameters is also provided. It is assumed that the 
communication system is asynchronous DS/CDMA 
MUD, over slow Rayleigh fading AWGN channel. 

The numerical results were obtained based on the 
averaging of 1000 simulation runs;and these results 
were attained in identical systems and channel 
conditions in order to provide fair comparison to other 
algorithms. Simulation parameters are as follow: 

TABLE 3 SIMULATION PARAMETERS 



Modulation 



BPSK 



Spreading code 



m-sequence 



Communication system 



Uplink asynchronous 
CDMA 



User number (K) 



20 



Channel 



AWGN with slow 
Rayleigh fading 



Path number (L) 



Path loss variance 



-5 dB 



The spreading sequences are selected as pseudo-noise 
(PN) m-sequence; in addition the number of active 
asynchronous users in the system is K=20; in all 
figures, it was assumed that the phases, amplitudes, 
channel gains and random delays of all users are 
perfectly known in the receiver, and users' power is: 



k ,1 



= 1, fork =\,2,...,K 



(16) 



In these simulations that will be seen, our proposed 
method is more effective and better in comparison to 
the two other algorithms, ACO and PSO. 

In Fig. 3, a comparison between the proposed 
algorithm with comprehensive search in the space of 
parameters known as ML versus BER of 
communication system is proposed. From this figure, 
it can be seen that increasing the length of m-sequence 
code causes the BER decreasing. Especially in code 
length equal to 10, the difference between our 
algorithm and ML is negligible. So using m-sequence 
with length 10 is suitable and it is used for the 
following simulations. Of course notice that codes 
with higher lengths have more complexity and are not 
suitable for our purpose of simple and low complex 
method. This simulation is done in SNR 15 dB. 




code length 

FIG. 3 THE BER COMPARISON BETWEEN ML AND PROPOSED 
ALGORITHM WITH DIFFERENT CODE LENGTHS (SNR=15 dB) 

The main parameter of each communication system, 
BER, is discussed in Fig. 4. This property of the 
proposed algorithm is compared to ACO and PSO 
algorithms and also to ML and the worst case as 
Matched filter. It is revealed that although in low SNR, 
there is no main difference among the methods; the 
proposed algorithm in high SNR converges to 
optimum MUD and the BER of this method is less 
than that in ACO and PSO methods. This figure shows 
the unsuitability of the matched filter detector. 




-d-MF 

-■-PSO 

-•-ACO 

— Sign_detector,code length=10 
-►-ML 



SNR(dB) 

FIG. 4 BER OF THE PROPOSED ALGORITHM IN COMPARISON 
TO OTHER METHODS VERSUS SNR. 

In Fig. 5, the effect of increasing the number of users 
on BER is analyzed for the classic algorithms and 
optimization-based techniques such as ACO and PSO. 
On the basis of this figure, it can be seen that the ACO 
and the proposed algorithm almost work similarly, 
and BER of these two algorithms does not 
considerably increase with increasing user numbers. It 
can also be observed that the performance of the 
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classic matched filter method is degraded substantially. 
This simulation is done for SNR 12 dB. 




9 10 11 

Number of users 

FIG. 5 BER OF THE PROPOSED ALGORITHM VERSUS NUMBER 
OF USERS IN COMPARISON TO OTHER METHODS. SNR IS 
FIXED AT 12 DB. 

In Fig. 6, the proposed algorithm with different code 
lengths is implemented and the numerical results are 
depicted. It is observed that with the code length 10, 
the algorithm converges to optimum MUD in high 
SNRs. Of course, for the other code lengths less than 
10 such as 8 and 6, the error is small and the 
convergence is acceptable to some extent. As seen 
from this figure, the proposed algorithm operates very 
well, and the convergence of this case is quite 
interesting. 




□ Sign detector, code length=6 
-•- Sign detector, code length=8 
-+- Sign_detector,code length=1 
-►-ML 



SNR(dB) 

FIG. 6 BER OF THE PROPOSED ALGORITHM VERSUS SNR FOR 
DIFFERENT CODE LENGTHS, COMPARED TO ML AND MF. 

Computational Complexity 

A common form in order to compare algorithms' 
complexities can be done through the O notation, 
which means the order of magnitude associated with 
the algorithm complexity. But comparing algorithms 



only with O can be insufficient. In order to express 
the complexity of the analyzed algorithms, it is 
essential to determine which instructions are carried 
out and how many times they are processed [1]. 
Suppose the number of transmitted bits is / and the 
number of users is K . For OMUD, the number of 
operations increases exponentially with the number of 
users, i.e. 2^ 7 . For our proposed algorithm, the code 
length and the size of correlation matrix of m- 
sequence codes are the most important factors of 
calculation. After that, the sign detector is very simple 
and has no complexity. The complexity of correlation 
matrix by code length 10 is fixed and it is calculated 
and stored before receiver in database of the system to 
use when needed. Fig. 7 shows that the complexity 
associated with our proposed algorithm is much less 
than ML and comparable to matched filter and ACO 
detector. 



^-MF 




-t-ACO 




—♦—Sign dete 


ctor,code lengths 


-►-ML 






20 

Number of Users 



FIG. 7 COMPUTATION ORDER OF THE PROPOSED METHOD 
VERSUS THE NUMBER OF USERS, COMPARISON TO ACO, ML 
AND MATCHED FILTER METHOD. 

Conclusions 

In this paper, multiuser detection based on a new 
algorithm was implemented, and through 
presentation of m-sequence codes and sign detector, 
the desirable method was achieved. When compared 
to sub-optimal algorithms such as ACO and PSO, the 
new algorithm shows better performance, and in some 
cases, approaches the optimal ML algorithm with 
minimum error. This algorithm is quite efficient, and 
shows acceptable performance in code length equal to 
10. Of course it should be noted that code lengths 
more than 10 has more complexity and more 
calculation time and codes with shorter length have 
not good performance as the 10-length one. 
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